Parameter estimation of binary compact objects with LISA: Effects of time-delay 
interferometry, Doppler modulation, and frequency evolution 

Aaron Roga and Sukanta Bosdil 
Department of Physics & Astronomy, Washington State University, 
1245 Webster, Pullman, WA 99164-2814, U.S.A. 

We study the limits on how accurately LISA will be able to estimate the parameters of low-mass 
compact binaries, comprising white dwarfs (WDs), neutron stars (NSs) or black holes (BHs), while 
battling the amplitude, frequency, and phase modulations of their signals. These modulations will 
arise from LISA's motion relative to the source, the evolution of the source emission frequency, and 
time-delay interferometry. We show that Doppler-phase modulation aids sky-position resolution in 
every direction, improving it especially for sources near the poles of the ecliptic coordinate system 
by a factor of at least six more than those near the ecliptic. However, the same modulation increases 
the frequency estimation error by a factor of over 1.5 at any sky position, and at a source frequency 
of 3 mHz, which is near the high-frequency end of the "unresolvable" white-dwarf confusion noise. 
Since accounting for Doppler-phase modulation is absolutely essential at all LISA frequencies and 
for all chirp masses in order to avoid a fractional loss of signal-to- noise ratio (SNR) of more than 
30%, LISA science will be simultaneously aided and limited by it. For a source with an initial 
instantaneous frequency z/ ^ 2.5 mHz, searching for its chirp or frequency evolution over a one- year 
duration worsens the error in the estimation of by a factor of over 3.5 relative to that of sources 
with ly ^ 1 mHz. Increasing the integration time to 2 years reduces this relative error factor to about 
2, which still adversely affects the resolvability of the galactic binary confusion noise. Thus, unless 
the mission lifetime is increased several folds, the only other recourse available for reducing the 
errors is to exclude the chirp parameter from ones search templates. Doing so improves the SNR- 
normalized parameter estimates. This works for the lightest binaries since their SNR itself does not 
suffer from that exclusion. However, for binaries involving a neutron star, a black hole, or both, the 
SNR and, therefore, the parameter estimation, can take a significant hit, thus, severely affecting the 
ability to resolve such members in LISA's confusion noise. Among the affected sources are galactic 
low-mass binaries containing a black hole, about which very little is known and which LISA may 
detect. Finally, we demonstrate how relative to a Michelson network a time-delay interferometric 
(TDI) data-combination network will have a different sensitivity sky-pattern and discuss its effect 
on parameter estimation. 

PACS numbers: 04.80.Nn, 95.55.Ym, 95.75.Pq, 07.05 Kf, 97.80.-d, 97.60.Jd 



I. INTRODUCTION 



The Laser Interferometric Space Antenna (LISA) is a proposed space-based broad-band gravitational wave (GW) 
detector funded by the National Aeronautics and Space Administration (NASA) and the European Space Agency 
(ESA) 01 . Unlike its extant and proposed earth-based counterparts, which include LIGO 2], VIRGO 3], TAMA 
H, GEO H, and LCGT 6], it enjoys guaranteed sources, specifically galactic binaries, such as AM CVn, HP Lib, 
WD 0957-666, RXJ1914+245, and 4U1820-30 In fact, LISA will be so sensitive that the population of galactic 

binaries of compact objects, which comprise white dwarfs (WDs), neutron stars (NSs), and black holes (BHs), will 
create source confusion noise that towers over LISA's instrumental noise near the sweet spot of its observation band 
1^,0,^^. This can potentially affect the detection of other GW sources. 

The confusion noise from WD binaries will be especially worse below about 3 mHz, where a good fraction of it, 
with low signal strengths, will be unresolvable. Galactic binaries that are bright enough or orbiting fast enough form 
a potentially resolvable population extending up to about 20 mHz. Our ability to resolve this latter class of sources 
will largely depend on how well we can estimate their signal parameters. Although the development of algorithms 
for resolving the component sources is still at its nascency [13, 113 , it is already clear that the smallness of the 
parameter errors is of utmost significance to the success of any of these algorithms. The focus of this paper is to study 
the role of different sources of these errors and determine the limits on how small we can keep them. 



* Electronic address: roganelli@wsu.edu| 
'''Electronic address: sukanta@wsu.edu 



2 



The problem of parameter-estimation accuracy has been studied for a variety of LISA sources in the past. Cutler 
did the first systematic study on this subject for monochromatic sources and super-massive black hole mergers 15j. 
Takahashi and Seto studied the problem exclusively in the context of low- mass compact binaries. They placed 
emphasis on the determination of the chirp parameter and its role in unraveling the distance to a binary, following an 
idea due to Schutz 17]. Barack and Cutler 18] studied this problem in the context of extreme-mass-ratio inspiraling 
binaries. Vecchio and Wickham 19] extended earlier parameter estimation studies to beyond the long-wavelength 
approximation by including the frequency-dependent modulation of LISA's beam-pattern functions arising at higher 
frequencies. They showed that invoking this approximation results in atmost a 3% loss in the SNR, but can have 
significant effects on the parameter accuracy (anywhere from 5% to a factor of 10) in the source frequency range 
3mHz <iy< lOmHz. 

All of these works, however, assumed that the data being filtered are two linearly independent LISA data streams, 
called Michelson variables, which measure the differential changes of length in LISA's three arms. While these early 
works were instrumental in deriving useful lower-limits on the parameter errors of high signal-to-noise ratio (SNR) 
sources, the attainability of these limits is unclear because the Michelson variables will suffer from the presence of 
excess noise arising from the laser frequency fluctuations and optical bench motion. The first noise source alone is 
sufficient to ruin LISA's characteristic strain sensitivity by several orders of magnitude. Fortunately, this critical 
problem has been addressed by devising a method for combining LISA's data streams in software to produce data 
combinations that have these noises mitigated to the level of LISA's noise budget This method is termed 

as time-delay interferometry (TDI) and will be introduced in Sec. ^ 

In this work, like Takahashi and Seto, we consider the parameter estimation accuracy of low-mass compact binaries 
that are not necessarily monochromatic. But unlike that work, we base our accuracy limits on TDI variables and 
do not use the long-wavelength approximation. We study the role of the chirp-parameter, which is tantamount to a 
monotonic increase in the inspiral frequency, in affecting the accuracies of all signal parameters. We also quantify the 
effect of the Doppler-phase modulation to the same end. Our study shows in what regions of the parameter space 
the effects of searching for the chirp parameter or the effect of Doppler-phase modulation are not negligible. This can 
be useful in simplifying prototypical studies for those sources for which these effects are small. On the other hand, 
parameter regions where a search for the chirp in a signal can hurt the accuracy of the other parameters estimates are 
identified. Moreover, we demonstrate that the sensitivity sky-pattern of a network of TDI variables is very different 
from that of a network of Michelson variables. A similar difference exists in parameter accuracies achievable by these 
two types of networks. Since any realistic search will be conducted in a network of the former type, inferences based 
on a network of the latter type will be off. 

The layout of the paper is as follows. In Sec. I^we briefly introduce a couple of TDI data combinations and show 
how they are obtained from combining the frequency shifts monitored in the laser beams traversing LISA's arms. 
Section ITIll describes the signal from low-mass compact binaries that are in circular orbits, but are spiraling in. The 
complete set of parameters specifying the GW signal from such sources is spelled out. A formalism for determining 
the limits on parameter accuracies in noisy measurements is given in Sec. IIVI This introduces the Fisher information 
matrix, whose inverse gives the limits on the elements of the variance-covariance matrix of parameter errors 23]. The 
effect of Doppler-phase modulation on the SNR of a signal is explained in Sec. and that of the source frequency 
evolution on the same quantity is explored in Sec. IVII The parameter errors themselves are studied as functions 
of the source frequency in Sec. IVIII and as functions of sky position in Sec. IVIIII We end with a discussion of the 
implications of this work in Sec. IIXI [2^ . 

II. THE DATA STREAMS 

The gravitational-wave information in LISA will be accessed through a total of twelve data streams exchanged 
among the three spacecraft, labeled i =1, 2, and 3 that are located at the vertices of an almost equilateral triangle, 
as shown in Fig. There, i and i* denote the two optical benches mounted in the ith craft. The arm-lengths of this 
triangle are labeled such that Li is the length of the arm facing vertex i. The unit vector, n^, specifies the orientation 
of the ith arm, and goes counterclockwise around the triangle in the figure. Six of these beams are inter-craft beams 
for monitoring gravitational waves, with 2 beams exchanged per spacecraft pair. The remaining six beams are intra- 
craft beams (not shown in the figure), with 2 streams exchanged between the two optical benches, i and i*, within 
each craft. 

An impinging GW will cause shifts in the central frequency of the laser beam exchanged between adjacent LISA 
craft. However, a GW is only one of the many sources of frequency fluctuations. There exist other noise sources that 
cause appreciable shifts in the laser frequency such that GW signals would be dwarfed in their presence. Shifts in 
the central frequency of lasers can be largely mitigated in hardware in a fixed arm- length interferometer. However, 
since LISA will be unable to maintain identical arm lengths while the craft formation fly around the sun, the beams 
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in each arm will possess significant frequency fluctuation noise. The second largest source of noise results from the 
mechanical vibrations of the optical benches on each craft. Time-delay inter ferometry has been shown to reduce both 
laser frequency and optical-bench motion noise to the level of LISA's noise budget, thus, allowing LISA to obtain the 
sensitivity necessary to function effectively as a gravitational wave detector 20] . 

To appreciate the basic idea behind TDI, we begin by introducing LISA's primary data products that will contain 
information about GW signals. Let /o be the central frequency of all the lasers in LISA, then the fractional shift in 
laser frequency in the beam originating at optical bench i can be written as 

Q{t)^^^ . (2.1) 

JO 

Likewise for the beam from the bench i*. It is important to note that what is directly measured is not the frequency 
fluctuation in each beam; rather it is the fractional frequency shift in a beam from the bench on one craft relative to 
that in the beam from the bench in an adjacent craft. This is how the six inter-craft data streams, labeled Ui and 
—Vi for "clockwise" and "counterclockwise" (see Fig. are obtained. The data stream Ui is obtained by beating 
the fractional frequency shift in the beam from bench {i — 1) relative to that in the beam originating at bench i* 
[25^ . For instance, Ui{t) = C3{t — I/2/c) — Ci{t). The remaining two streams, namely, U2 and Us^ can be obtained by 
cyclic permutation of the indices in the Ui expression. In the rest of the paper, we set the speed of light, c, to unity. 
The —Vi data streams are obtained in a similar fashion by measuring the fractional frequency shift in the beam from 
bench {i + 1)* relative to that originating at bench i. The time-shift operator d is defined by its action on a data 
stream x{t) as p26j: 

Cix{t) = x{t - Li) , (2.2) 

where the label i denotes the arm along which the time-shift is affected. Then, Ui{t) can be written alternatively as 
Ui{t) = C2Cs{t) — Ci{t) and, similarly, for the other streams. 

The intra-craft beams exchanged between adjacent optical benches on every craft provide additional information 
about optical-bench motion noise that is common to multiple data streams. Although there are a total of 6 intra-craft 
beams, the two beams per craft can be combined in such a manner so as to produce a single data stream containing all 
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the relevant information regarding these noise sources. This can be seen easily by examining a single craft's internal 
data streams. At craft 1, by beating the fractional frequency shift in the beam from bench 1 relative to that in the 
beam from bench 1*, one forms the data stream, 

Wi = {Ci - Us • vi) - (Ci* + n2 • vi*) + na • ui + n2 • ui* . (2.3) 

where Ui^i* are the random velocities of the proof masses on benches i and z*, respectively, and Vi^i* are the random 
velocities of the optical benches themselves. Two other intra-craft data combinations, W2 and W3, can be obtained 
by the cyclic permutation of indices in the above expression. 

There exists a set of polynomial operators, pf^ qf, and r^, that when acting on the nine basic data streams (viz., 
six inter-craft streams, Ui, Vi, and three intra-craft streams, Wi) allow one to construct several time-delayed data 
combinations that are free from the dominating noise of the laser-frequency fluctuation and optical bench motion. 
This technique first developed by Tinto and Armstrong 20j is referred to as time-delay interferometry (TDI) and the 
resulting data streams are known as TDI variables. Subsequently, it was shown by Dhurandhar et al. 27] that a 
host of such laser-frequency canceling streams exist and form the module of syzygies. In an earlier work 26], it was 
shown that any of these data combinations, or pseudo-detectors, obtained from operating with these polynomials on 
the basic streams can be expressed as: 



Trace [e^^ • Z] , (2.4) 



where A is the pseudo-detector index and 

-^=llt- qi q'i 1 and Z = { V2 U2 W2 ] . (2.5) 
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Although there are many sets of noise-canceling pseudo-detectors, one well-studied set of pseudo-detectors is defined 
by the choice H,!!!^]: 

/ 1 Cs C2C3 \ / C1C2 1 Ci \ 

= 1 C1C2 C2 , = Cs 1 C2C3 , (2.6) 

V1 + C1C2C3 C3 + C1C2 C2 + C1C1/ VCs + CiCs I + C1C2C3 C1+C2C3/ 

and 

/ C2 C2C3 1 \ 

= C1C3 Ci 1 • (2-7) 
VC2 + C1C3 C1 + C2C3 I + C1C2C3/ 

These pseudo-detectors were termed as a, and 7, respectively, by Armstrong et al. [2lL[2^. 

It is important to note that the data combinations above are not statistically independent. However, it is possible 
to create noise- independent pseudo-detectors by diagonahzing the noise-covariance matrix of the above combinations. 

A=^{-a + j) E = ^{a-2/3 + j) T = ^{a + p + j) (2.8) 



One such triplet is [2 



where T is akin to a Sagnac response, which is devoid of a GW signal. The data combination A is not to be confused 
with the pseudo-detector index A: Unlike the former, the latter will appear only as a superscript or subscript of 
another symbol. What is most important to note is that the eigenvalues of the above-mentioned noise-covariance 
matrix are degenerate, with A and E both having the same noise variance. Therefore, any pair of pseudo-detectors 
that lie in the plane normal to T and are orthogonal to each other will be noise-independent as well. An infinity of 
such pairs are obtainable by simple rotations about T. Unless specified otherwise, for the remainder of this paper, 
the pseudo-detector triplet or network considered is 

A = cos(7r/3)A - sin(7r/3)£; = x^ E = - sin(7r/3)A - cos(7r/3)£; = x^ f = T = x^ . (2.9) 

An alternative data combination that has been studied in the literature is the pair of Michelson variables 

h'{t) = {e,{t) - e2it)) , 

h"{t) = {ei{t) + e2{t) - 2e3{t)) , (2.10) 
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where Si is the strain in the ith arm. These variables suffer from the frequency fluctuation and optical bench motion 
noise and will not be useful in actually detecting signals once LISA is operational, except at very low frequencies. 
However, owing to their simple relation to GW polarization components 0, they have been found to be a useful 
starting point for data analysis studies. 

As discussed in Ref. |2^J , the three noise-independent flrst generation TDI variables can be expressed as 

x^{t)=n^{t)^h^{t) , (2.11) 

where 

h^it) = E [Pfv^^it) + QfU^^it)] (2-12) 

i=l 

is the gravitational- wave signal in the Ath pseudo-detector and n^{t) is the remaining instrumental noise in it. (Note 
that the sum on the right-hand side adds up to zero for A = 3.) The noise is assumed to have a Gaussian probability 
distribution with a zero mean. The noise covariance matrix elements are given as follows: 

n^*(/)n^(/') = \p^^Hmf - f)S^^ , (2.13) 

where P^^\f) is the one-sided noise power-spectral density (PSD) of the Ath pseudo-detector. 
The noise PSDs for the above mentioned pseudo-detectors can be expressed as Iscj : 

p(i)(j) = p(2)(/) = 8sin2(7r/L){2[3 + 2cos(27r/L) + cos(47r/L)]PP^°°^ 

+ [2 + cos(2^/L)]P^h°*} , 
P^^\f) = (2 + 4cos'(2^/L)) [(4sin2(^/L)pP^°°^ + P^^°*] , (2.14) 

where the noise PSDs arising from the photon-shot noise and the proof-mass noise are estimated to be P^^°* = 
1.8 X 10-3^[//1Hz]2Hz-^ and PP^^^^ = 2.5 x 10-'^^[//1Hz]-2Hz-\ respectively Q. The proof-mass noise affects 
LISA'S sensitivity at low frequencies, while the shot noise affects it at higher frequencies. Note that P^^\f) is the 
noise PSD for pseudo-detectors A as well A; similarly, P^'^\f) is the noise PSD for pseudo-detectors E as well and 
so on. 

The TDIs introduced above are the so-called flrst generation TDIs. These TDIs assume that the light travel time 
for data stream Vi and data stream are identical. Realistically, however, the rotational and orbital motion of 

LISA would prevent the noise contribution of the laser-frequency fluctuations from being mitigated to the level of 
the secondary noises |^ Is^, • In order to tackle this problem, new pseudo-detectors were introduced as simple 
differences of their flrst generation counterparts, appropriately time-shifted: [3^ : 

h^{t)=h'^{t)-h'^{t-3L) , (2.15) 

where the subscript 2 indicates that these are second generation data combinations. One can also deflne the second 
generation noise PSD as: 

pa),(2),(3)(^) = 4sin2(27r/L)p(i)>(2).(3)(/). (2.16) 

Again, the subscript 2 denotes that these are second generation noise PSD's. The flrst and second generation noise 
PSDs are compared in Fig. [21 

In the rest of the paper, we will discuss parameter estimation of binary inspiral signals based on the flrst generation 
data sets (unless otherwise specifled). The above transformation can be implemented straightforwardly to get the 
estimates for the second generation data sets. 



III. THE GRAVITATIONAL- WAVE SIGNAL 



Consider a GW source located at a distance R relative to the solar system barycenter (SSB). Let (O^cj)) denote 
the sky position of the source in the ecliptic coordinate system centered at the SSB. Then the sky-position vector is 
given by w = (sin 6> cos 0, sin 6> sin ^, cos 6>). In the transverse-traceless gauge j35(J, the space-space part of the metric 
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FIG. 2: Plots of the first and second generation noise power-spectral densities (in units of Hz ^) 



perturbation due to a gravitational wave emitted by this source and observed at the spacetime location (t, r) can be 
expressed as 



hmn{t, r) = - W • r){OmOn 



■hx{t-w - r)(6>^ 



(3.1) 



where m and n are spatial indices, Om and (j)n are components of unit-normal vectors and 0, respectively, that 
define a transverse plane with respect to the wave-propagation vector, — w, and /i+ and hx are the two polarization 
components of the impinging wave. Therefore, the strain in LISA's ith arm is characterized by the contraction of the 
arm's "detector" tensor, n^n^^ with the perturbation matrix hmn{t,'r): 



where 



h,(t) = h+m+{e, <p) + h^ (t)eix {e, 0) , 

^,+ {t) = {e-Ui{t)f -{^■Ui{t)f and CixW =2(^-n,(t))(0-ni(t)) , 



(3.2) 
(3.3) 



define that arm's beam pattern. 

In this paper, we consider how well one would be able to estimate the parameters of a nearly monochromatic source, 
such as a low-mass compact binary, in the LISA data. Let the masses of the two stellar members of these sources be 
mi and m2. Then, in the quadrupole approximation, the amplitude of its signal at the SSB can be expressed as 



H{uj) = L188 X 10" 
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(3.4) 



where M = rr^-l^rri^^ l{m\ + m2)^l^ is its chirp mass, R is its luminosity distance, and (jJ = 2i^v is its angular 
frequency. 

For a sufficiently high-mass binary, the inspiral of its members will cause the source frequency to increase perceptibly 
in the LISA band. On the other hand, for small chirp mass and emission frequency, as studied here, the rate of 
frequency evolution is approximately ^] 



48 (GM\ 
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(3.5) 
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where G is the universal gravitational constant. 

Both the amplitude and the phase of a GW signal will undergo modulation due to LISA's changing orientation and 
motion relative to the source. The amplitude modulation is determined by the sky position angles and the source 
orientation. The latter can be specified in terms of the polarization angle, ?/^, and the inclination, ^, of the binary's 
orbit relative to the line of sight. As was shown in Ref. 26], this modulation in the AXh pseudo-detector, arising from 
LISA'S zth arm, is captured by the complex extended beam pattern functions, 

E,^:=r/D,f , (3.6) 

where the sum over the values of p = ±2 is implicit. The advantage of the above decomposition of the extended beam 
patterns is that all the source-orientation dependence is limited to the Gel'Fand functions [33. 13^ 13^. 

T2^^ (V^, ^, 0) = ^ (1 ± cos if exp (Ti2V^) • (3.7) 
The infiuence of the source location and the arm orientation on the amplitude modulation is factored in 

^^2i = ^|M^lte+TiC,x) , (3.8) 



M- = ,sinc Lu^l±^) - .sine U^^^^") . (3.9) 



where 

As was shown in Ref. the strain in the Aih. pseudo-detector can be expressed as: 
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h^{t)=H{u,)Y,^[e-''Et*{t)St{t)\ , (3.10) 

i=l 

where 5 is the initial phase and 

= ^iu;it-Liri)+ii^d;t^+af)+i(l)Dit) (3.11) 

describes the phase evolution of the signal. Above, af = arg(M/^) denotes the modulation in the phase arising from 
LISA'S changing beam patterns, and 

o ' (3-12) 



2V 

where L^w • ri/V^ is the delay in the arrival of the signal at the ith craft relative to LISA's centroid. Finally, 

<^i:)(^) = ^<^0 sin^cos [ — ^ J (3.13) 



gives rise to the Doppler-phase modulation owing to LISA's motion relative to the source. Here, (= 1 AU) is the 
average distance to the sun from LISA's centroid and P (=1 year) is the time taken for LISA to complete one full 
orbit around the sun. 

The strain expression given in Eq. (|3.10|) will be used below to ascertain how well we can estimate a slightly 
chirping {co ^ 0) binary's parameters with LISA. The above discussion shows that such as signal will be characterized 
by eight parameters, 'd = {T^z>, i^, ^, V^, cos ^, ^, cos 6>, Ti^}, such that, "d^ = T^z>, i^^ = H, and i^^ = Tu. Note that 
all parameters are dimensionless. 



IV. ESTIMATING SIGNAL PARAMETERS 



In estimating the parameters of a signal, an observer has to confront the inherently noisy nature of the detector 
or receiver, which introduces a degree of randomness in the parameter measurements. Thus, an ensemble of multiple 
copies of the detector will return a distribution of parameter values that is affected by the characteristics of the 
detector noise. The spread in the distribution indicates the noise-limited accuracy of the detector. 
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Before we define our estimator and determine its accuracy, note that even in the absence of instrumental noise, a 
typical detector has physical limits on its measurement accuracy. In this context, it is important to note that there 
are two important time-scales that affect LISA's observations: The first arises from LISA's orbital baseline, the mean 
light-travel time for which is close to 1000s. Thus, the frequency of a wave with wavelength equal to this distance is 
about ImHz. The other scale arises from LISA's 5 million-km arm-length, which corresponds to a light-travel time of 
16.67s. Thus, waves with frequencies that are odd-integral multiples of 0.03Hz will drive the two craft at the ends of 
an arm out of phase by 180 degrees. This makes LISA's sensitivity oscillatory at the higher end of its band (see Eqs. 
I^TT^ and Fig. □) 

To ascertain how large the noise-limited errors will be in the parameter values, we take the parameter values 
themselves to be the maximum likelihood estimates (MLEs). Owing to noise, the MLEs, i^, can deviate from the true 
values of the signal parameters and can influence the determination of other parameters. The magnitude of these 
deviations and influences can be quantified by the elements of the variance-covariance matrix, = {'d^^'d^} j2^ . 

A relation between the and the signal is available through the Cramer-Rao inequality, which dictates that 

II 7 II > lirir^ , (4.1) 

where T is the Fisher information matrix: 

2 



2 7^ 

^ E pT^) I [dmh^it-^ ^)] [9nh^(t-^ ^)] dt ■ (4.2) 

In this paper, we determine the values of Tmn for signals from low-mass compact binaries in LISA data. Therefore, 
Frnm givcs the lower bound on A'd^ (which is the expected random error in the MLE of i^^). The two are equal in 
the limit of large SNR. 

Since the amplitude and the initial phase are extrinsic parameters, one can analytically maximize the likelihood 
ratio with respect to them |2^. Here, we assume that this has been done. In such a case, the Fisher information 
matrix is six-dimensional. The errors in the phase and the amplitude can be derived in terms of the errors in the 
remaining six parameters as demonstrated, e.g., in Ref. This is the method we use to compute the error in H. 

One can similarly infer the error in 5. 

The errors in the sky-position angles will be presented in terms of the error in the measurement of the sky-position 
solid angle, defined as: 

AQs = 27r (A^^A^^ - (A^^A^^)) . (4.3) 
We define the error in the measurement of the source-orientation solid angle similarly: 

AQl = 27r {A^^A^^ - (A^^A^^)) , (4.4) 
where the subscript L denotes the orbital angular momentum of the binary. 



A. Scaling of parameter errors with SNR 

Although our expressions above can be used to compute the error estimates for any parameter values, we illustrate 
our results for specific values of some of these parameters: Unless specified otherwise, we use an integration time of 
one year, an emission frequency of u = 3 mHz, and source orientation angles ip = tt/S and t = 7r/4. The emission 
frequency is chosen to be large enough so that the Doppler-phase modulations have a significant effect. The source 
orientation is chosen somewhat arbitrarily, except that we ensure that the signal is not linearly polarized. 

All the parameter estimates studied below are obtained by normalizing the estimates given in Eq. (|4.2|) by the 

r 2 4 4 1 "^ 

network sensitivity, , ^ • Thus, the parameter-estimate values shown in the plots here are for a 

signal with an SNR of 1. To get an estimate for an arbitrary SNR, one simply multiplies the value shown in the plots 
by the SNR scale given in Table U 

An alternative way of understanding the parameter estimate normalization is through the Fisher information matrix. 
The idea here is to first divide the Fisher information matrix by the network template norm, Y1'a=i{^'^ ^ h"^){A)^ 
follows: 

(4.5) 
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TABLE I: The table lists how the errors in different parameters scale with the SNR. Since all our error plots are given for an 
SNR of unity, to assess the error for any other SNR one simply needs to multiply the plot value with the scale read from this 
table. 
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One then obtains the normalized error estimates from the covariance matrix derived by inverting F^^- 

V. DOPPLER-PHASE MODULATION 

When a template does not track the Doppler-phase modulation of a signal, over time it gets increasingly phase- 
incoherent relative to the latter. We quantify the degree of coherence between the two by the fitting factor: 



(5.1) 



where the value of the "template parameter" i^' is, in general, different from that of the "signal parameter" t^. The 
plot of the above function, with h^{'d') as the template without Doppler shifting and h^{'d) as the signal with that 
shifting accounted for, is given in Fig. |31 As can be inferred from that figure, when the Doppler-phase modulation 
is not accounted for in the search templates, the SNR suffers an appreciable drop across the entire LISA band. It 
also exhibits large oscillations with a period that depends on LISA's orbital speed and that decreases linearly with 
increasing source frequency, as is characteristic of Doppler frequency shifting. It is important to note the locations of 
the nodes in the plot: The first minimum is at approximately / =1/2L, which occurs when the phases of a gravitational 
wave impinging along an arm are different by 180° at the two craft at the two ends of that arm. 




FIG. 3: This is a plot of the fitting factor between two similar templates, one with the Doppler-phase modulation included and 
the other without it. 



In conclusion, since the effect of discounting Doppler-phase modulation is so marked on the SNR at any LISA 
frequency of interest, one must account for its effect on parameter estimation throughout the LISA band. 
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VI. EFFECT OF SOURCE FREQUENCY EVOLUTION 

For binary compact object that are spiraling in fast enough, LISA wih be able to measure their chirp mass and 
luminosity distance. Expression (|3.4|) shows that the amplitude, H{u;)^ of a signal from such an object depends jointly 
on r and Aic- However, neither of these parameters affects any other part of a monochromatic (or z>) signal, as given 
in Eq. (|3.10|) . Thus, one can not separately estimate either of them purely from the amplitude. Nevertheless, if a 
signal has an appreciable amount of inspiral or frequency evolution, defined in Eq. (|3.5j) . then the measurement of z>, 
along with H{uj), determines both the chirp mass and the source distance. 

In order to ascertain the individual masses of the binary, additional information is necessary. Traditional astronom- 
ical methods using optical or radio measurements can be used to help identify the total mass of the binary system. 
However for many unknown sources or optically invisible sources this is not an option. In some cases where the binary 
orbit is noticeably eccentric, the accurate tracking of the emission frequency can provide information on the total 
mass of the system. This datum, together with the chirp mass, can be used to infer the individual masses of the 
system [i^ . 



1.2 




FIG. 4: Fitting factors of low-mass compact binaries with a range of chirp masses are plotted as functions of frequency. Here 
the fitting factor is between two template families, one with a physical set of eight parameter values and another with the 
chirp set to zero, artificially, for all source frequencies. For the above plots, the sources are located on the ecliptic, with 
= 7v/2, (/) = 7r/2, and their orientation is ^ = 7r/4, -0 = 7r/3. As a reference, an equal-mass binary, with mi = m2 = 0.35 M© 
has A4c — 0.3 M©. Fitting factors for equal-mass binaries with m =1,2, 3,... 6 Mq can also be read from here. 

Before analyzing the quantitative effects of a chirping waveform on parameter estimation we briefly study the role 
that a non-negligible z> plays in searching for a signal. Speciflcally, we use the fltting factor deflned in Eq. (|5.1|) 
to determine the fractional loss in the SNR while searching for a chirping source with monochromatic (i.e., z> = 0) 
templates. When subtracted from unity, the above function measures the fractional drop in the SNR owing to a 
parameter mismatch of Ai^ = 'd — 'd' . In Fig. ^ we plot m{'d^'d') with /S.'d'^ = 0, for all m 7^ 0; the "template 
parameter", 1^^^, is set to zero, while the "signal parameter", i?^, is determined by Eq. (|3.5|) for any given emission 
frequency and chirp mass. For any given chirp mass, note how rapidly the SNR drops as a function of frequency. 
Also, since at any given frequency the chirping causes a smaller mismatch between the two template families for a 
smaller chirp mass, the SNR drop occurs at a higher frequency for a smaller chirp mass. Just like Fig. |31 this fltting 
factor too has oscillatory features. However, owing to the polynomial frequency behavior of z>, the oscillations here 
do not exhibit the linear decrease in the node-spacing found in Fig. |31 

We now turn our attention to the role of z> in the estimation of the signal parameters. As shown in Sec. II VI the 
expansion of the template parameter space in order to include another parameter affects the parameter estimation 
through the Fisher information matrix, which is a different construct than the fltting factor. Once a detection has 
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been made, the error values presented here provide the standard deviation that should be expected for an ensemble of 
signal measurements made by LISA. However, one still needs to establish for what ranges of co and Aic it is possible 
to measure a non-vanishing chirp in a signal. Outside these ranges, including u in ones template-parameter space 
can result in worse estimates for other parameters. A smaller parameter volume for estimating parameters is also 
desirable because it leads to a lesser computational burden: This is more so since the problem of detection in the 
presence of LISA's source confusion noise is intricately linked with the problem of parameter estimation. 

10^ 



10^ 
10" 



*t: 10"' 
10-^ 



'"id"" 10"' 10"' 

f(Hz) 

FIG. 5: These plots illustrate the values (in number of cycles) of the chirp parameter for two values of the chirp mass 
(y\4c = O.3M0 and S.SM©), along with the errors in determining the chirp and frequency for an SNR of 10. All plots are for an 
observation period of T = 1 year. Note how for the larger chirp mass the error in determining the chirp remains comparable 
to the value of the chirp itself up to about 1.5mHz, before improving at higher frequencies. 

Figure [51 shows that for binaries with Aic — 5.3 M©, the error in the estimation of its chirp is at least as large as 
the value of the chirp itself (for an SNR of 10) for u ^ 2.5 mHz. As shown in Fig. |H1 this is very close to the frequency 
where the inclusion of the chirp starts hurting the estimation of the frequency. On its own merit, this would argue for 
the dropping of the z> parameter from the search templates (i.e., setting = 0) for Mc ^ 5.3 Mq and u ^ 2.5 mHz. 
However, the fitting factor plot shows that doing so can vastly reduce the probability for detecting such signals with 
v ^ 1 mHz. We have thus, identified a region in the parameter space where entertaining the prospects of a detection 
requires taking a beating in estimating the source parameters. This adds an additional layer of challenge in the quest 
for cleaning what is deemed as the resolvable part of the confusion noise spectrum. 

VII. FREQUENCY DEPENDENCE OF PARAMETER ERRORS 

Imagine a one-parameter family of sources, which have everything else identical, except for the emission frequency, 
u. The parameter-estimation errors for such a family are not invariant across LISA's bandwidth. This is because both 
the Doppler-frequency shifting and the frequency evolution increase with u. Additionally, the sensitivity of LISA itself 
varies across the band, as shown in Fig. [7\ For nearly monochromatic signals, the latter effect gets essentially factored 
away when considering S NR-normalized estimdites. Figure [SI shows the frequency dependence of the parameter errors, 
normalized for an SNR of unity, for four different cases: Case (a) depicts these errors when a signal has no Doppler- 
phase modulation, and no frequency evolution (where the latter is unphysical). Case (b) presents the same when 
there is frequency evolution natural to the binary, but no Doppler-phase modulation (such as if LISA were at the 
SSB). Case (c) includes the effect of Doppler phase modulation but has the frequency evolution (artificially) turned 
off. Finally, plot (d) has effects from both the Doppler phase modulation and the frequency evolution. 

It is interesting to note that in Fig. the error in the estimate of each of the parameters is essentially independent 
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[a] [b] 




FIG. 6: These figures show how the errors in the different signal parameters behave as a function of the source frequency 
for the following four cases (a) with no Doppler-phase modulation and no chirp, (b) with chirp only, but no Doppler-phase 
modulation, (c) with Doppler-phase modulation, but no chirp, and (d) with both Doppler-phase modulation and chirp included. 
The sky-position assumed is: = 7r/2, = 7r/2; also, ^ = 7r/4, V^^tt/S. 

of the emission frequency. The primary contribution to the frequency dependence of the Fisher information matrix 
in Eq. (|4.2|) arises from LISA's own sensitivity reflected in its noise PSD. The elements of this matrix depend on 
P^{f)i which is the same for both A = 1,2. However, when the parameter errors are obtained by scaling this matrix 
by the template norm, as described in Eq. (|4.5|) . the factor of P^{f) for near- monochromatic signals gets canceled 
away, thus rendering the plots in Fig. El mostly flat. 

We next consider the effect of Doppler-phase modulation on the parameter estimation accuracy. Comparing Figs. 
andEfc shows that the errors in all parameters are mostly unaffected except for the error in sky-position. That error 
actually worsens with frequency in the 0.1-1 mHz band, before rapidly improving at higher frequencies. This is because 
at wavelengths longer than 1 AU, the Doppler-phase modulation acts like noise on the sky-position dependent cr^ term 
in the phase of the signal in Eq. (|3.11|) . thus, affecting sky-position estimation. On the other hand, at high frequencies, 
it does a better job at resolving a binary's direction than do the beam-pattern-function induced modulations of the 
amplitude and the phase. 

The effect of the chirp can be understood by comparing Figs. Efc and Eli (or, alternatively. Figs. Eh andEb): 
It worsens the estimation of the amplitude and source orientation at low frequencies, i.e., for ^ 1.5 mHz. It also 
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FIG. 7: These are plots for the sensitivity of LISA as a function of the frequency (in Hz) for the first and second generation 
TDIs, respectively. 



reduces the accuracy of its own estimate and that of the source frequency at high frequencies, i.e., for u ^ 1.5 mHz. 
Here, one may wonder why at frequencies below 1.5 mHz, the chirp's inclusion ruins the determination of the source 
orientation, but not the frequency. The answer is that during the course of a year one integrates through a very 
large number of cycles of a signal, which averages out the noise in the estimation of the frequency arising from the 
chirp. The same would be true about the orientation if one were to integrate for a large number of LISA orbits (or 
years). Interestingly, at frequencies higher than 2 mHz, while the chirp leaves the estimation of the source orientation 
unaffected, it hurts the estimation of the source frequency. This is because starting at around a milli-Hertz, the chirp 
acts as a source of noise in the determination of the frequency. 

In summary, the general behavior of the SNR normalized parameter-estimation errors as functions of u can be 
understood as follows. Primarily, there are two competing parts of the signal that provide information about the 
source: The amplitude modulation and the phase modulation. The amplitude is modulated differently throughout 
the course of a year depending on the source location and orientation. The amplitude modulations provide information 
about these two source parameters in a manner that the SNR-normalized error in them remains invariant across LISA's 
band (as can be inferred from Fig. El a). Whereas tracking the phase provides information about the source location, 
source frequency, frequency evolution and, hence, the chirp mass. The main difference between the two modulations is 
that the phase modulation is highly dependent on the emission frequency of the source, whereas amplitude modulation 
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is not. Since the amplitude modulation is less sensitive to the emission frequency, its contribution to any parameter 
estimate is relatively unchanged over LISA's band. Contrastingly, phase modulation provides less information at 
low emission frequencies. This is due to the relatively fewer cycles one obtains in a given observation time for low 
frequencies compared to higher ones. For the frequency range that LISA is projected to observe, this translates into 
roughly a factor of 10^ less cycles that would be obtained for the smallest observable frequencies with respect to the 
largest ones. However, as the frequency increases the contribution from the phase modulation eventually surpasses 
that due to the amplitude modulations. 

A. Excluding frequency evolution from the templates 

The preceding discussion shows that at any given LISA frequency searching for the frequency evolution of a binary 
affects the determination of one or the other parameter, except the sky position. For a source with an initial frequency, 

^ 1.5 mHz, its inspiral hurts the error in the estimation of z^ by a factor of over 3.5 relative to that of sources with 
V ^1 mHz. We have verified that increasing the integration time to 2 years reduces this relative error factor to 
about 2, confirming what is shown in Ref. [la|- This reduction may be insufficient for the resolution of the galactic 
binary confusion noise. This makes a case for a mission lifetime that is longer than 2 years. Otherwise, the only other 
option available for maintaining the parameter accuracies is to exclude the chirp parameter from the search templates, 
since doing so improves the SNR-normalized parameter estimates. This, however, has different consequences for the 
parameter accuracies of binaries with different chirp masses. 

The fitting-factor plots in Fig. 2|show that for the lowest astrophysically interesting chirp mass, Mc — O.SM©, such 
an exclusion has little effect on the SNR until about 4 mHz, which happens to be just below the high-frequency-end 
of the confusion noise arising from extreme mass-ratio inspiral (EMRI) sources. The effect on the estimates of the 
other parameters is a systematic bias by an amount that depends linearly on T^z> 43], which is about 0.1 (cycles) 
for T = lyear and at = 4 mHz. While this introduces a negligible bias in the estimates of all parameters, it has 
the advantage of reducing the level of the random estimation errors to the same level as that shown in Fig. Et- At 
frequencies as high as even 6 mHz, the SNR drop due to such an exclusion would be more than 30%, thus, jeopardizing 
both the prospects of detection and the adequate parameter resolution required for cleaning the WD- WD confusion 
noise. 

On the other hand, note that for a binary with Mc ^ 5.3M0, such an exclusion causes the same fractional drop 
in the SNR at frequencies as low as 1.5 mHz, which is where the parameter estimates start getting affected (see Fig. 
int> or|nii). Thus, for binaries with Mc > 5.3M0, turning off the chirp parameter in the search templates is unlikely 
to improve the parameter estimation for u > 1.5 mHz. Among the sources most affected from this "chirp search" 
problem are galactic low-mass binaries containing a black hole. 

VIII. SKY-POSITION DEPENDENCE OF PARAMETER ERRORS 

In this section, we study the behavior of the parameter errors as functions of the direction to a binary. We illustrate 
the sky-position dependence of the SNR-normalized parameter errors for the same source polarization state and 
frequency as the one chosen in Sec. IVIII for the frequency plots in Fig. El The source frequency choice, viz., u = 3 
mHz is interesting for multiple reasons: As observed above, it is close to the high-frequency-end of the unresolvable 
WD- WD confusion noise. Figure shows that it is also in the region of transition in the behavior of the parameter 
errors as functions of frequency. All the same, it is not too small for the Doppler-phase modulation to have a negligible 
effect on the parameter accuracies. 

When studying these errors as functions of frequency in Fig. [SI we had fixed the sky position of the source. For 
this reason, while computing the SNR-normalized errors, we did not have to worry about the possibility that the SNR 
may vary with the sky position. We begin by addressing that possibility below. The dependence of the SNR on the 
polarization state can be studied similarly. 

A. Sky-position dependence of the SNR 

The SNR normalization performed on the Fisher information matrix in Eq. 14.51 is useful in that the parameter 
errors computed from such a matrix can be easily scaled for a signal of any SNR using the scaling factors given in 
Tabled In this context, it is important to note that the SNR of a standard candle is not uniform across the sky owing 
to the non-uniform sensitivity of LISA to different sky positions. We plot this quantity for four different cases in Figs. 
IHlandini where the latter set is the 3D rendition of the former set. These plots, called "sky plots", are normalized so 
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FIG. 8: These are plots of the template-norm-squared for the same four cases of LISA motion and source frequency evolution 
considered in Fig. El The plot for each case is normalized to have a maximum value of unity. Both and (j) are given in 
radians, with the origin, (0, 0), in the bottom left corner of each plot. Note that this is the first of a series of colored sky-plots; 
gresy scale prints of these plots can give a misleading impression of the color values. 



that the maximum value attained in each case is unity. These four cases are identical to those considered for studying 
the frequency dependence of parameter errors in Figs. As in that study, here too Fig. [Ht (or, equivalently. Fig. 
Et) is for an artificial signal, with z> set to zero), when compared with Fig. [Hb, it allows one to graphically assess 
how useful or detrimental the presence of frequency evolution is in detecting a signal. Similarly for Fig. |Hfc, when 
compared with Fig. [Hli. Note that as per the definitions of cases (a) and (b). Figs. [Ht andlHb have the Doppler-phase 
modulation turned off. 
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FIG. 9: Spherical renditions of the template- norm- squared plots shown in Fig. |Sl 



B. Explaining the sky-position dependence 

Given that LISA's orbital motion is confined to a nearly circular orbit on the ecliptic, one may expect the SNR 
plots to be axisymmetric. This is indeed the case for the network comprising the two Michelson variables, as can be 
seen in Fig. El where we limit our focus to their geometric sensitivity and momentarily ignore the fact that they will 
be severely limited by the laser frequency- fluctuation noise. Figures ITUk andlTHb reveal that although the SNR for 
neither of the Michelson variables is axisymmetric, their quadrature combination is so. For comparison, we present 
the analogous plots for the A and E TDI variables in Figs. We normalize all these plots relative to the maximum 
value attained among the template norms of the individual Michelson and A-E variables. 

Interestingly, however, the SNRs for the A and E TDI variables, plotted in Fig. ^2 are more dumbbell-shaped than 
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axisymmetric, where the pattern depends on the initial orientation of LISA. Also, since there is a slight difference 
between the overall sensitivities of the two variables, their quadrature combination is more of an inflated version of 
the A variable. As one would expect, it also resembles well the SNR of Fig. [HK- 

Figure El illustrates the sensitivity of a single arm to all sky positions. What is evident from this figure is its 
striking similarity to the A TDI data combination. The features of this plot are well understood geometrically. For 
a single arm, spinning and orbiting the SSB, there exists a periodicity of tt radians under axial rotations. However, 
this figure shows that there are two distinct and orthogonal directions, and their antipodes, that stand out. These 
locations correspond to the global maxima and minima in the sensitivity of this arm to impinging gravitational waves. 
These very same extrema are present in the SNR plots as well. To understand the origin of these global extrema, the 
orbital motion of LISA's centroid is not as important as the spin of LISA's individual arms. Therefore, consider how 
a single arm spins during the course of a year. Let it begin in an orientation that is parallel to the ecliptic plane. 
In this orientation, a low-frequency source that is in a direction perpendicular to the arm will create the maximum 
strain in that arm. Conversely, when such a source is in a direction along the arm, the strain caused in that arm 
is zero. Let us call these directions 1 and 2, respectively. With the exception of their antipodes, there is no other 
orbital position that will attain such extrema in the sensitivity for that particular arm during its orbit. As the arm 
rotates out of the ecliptic plane, the sensitivity to direction 1 will begin to degrade. However, because direction 1 
will never be parallel to the same arm, it will never attain the same minima that direction 2 attained initially and, 
vice versa. Therefore, the orbital position where these global extrema are obtained is forever imprinted upon the 
integrated results regardless of the duration of observation. 




FIG. 10: The template- norm- squared for the Michelson variables as functions of sky position: for Michelson variable I as a 
function of sky position (left), II. for Michelson variable II as a function of sky position (center), III. for the network comprising 
both Michelson variables (right). All plots are presented for the following source parameters: {^ = 7r/4,'0 = 7r/3,z^ = 3 mHz} 



C. Parameter errors 

We plot the parameter errors as functions of sky position in Figs. ^1^^ and IT^ for the same source polarization 
state and frequency as described earlier in this section. Note that these parameter-error sky plots are presented for 
the AET pseudo-detector network. Altering the polarization state, or the pseudo-detector network will, in general, 
change the details found in these plots. However, certain qualitative aspects, discussed below, will remain unaffected. 

While all four cases (i.e., (a)-(d)) are plotted for sky-position and frequency errors in Figs. H^andlT^ respectively, 
only the full case (i.e., case (d)) is plotted for the remaining parameters in Fig. This is because the plot for each 
case of these remaining parameters strongly resembles the frequency-error sky plot of the corresponding case in Fig. 
fm Indeed, a quick examination of all the error plots reveals that they can be divided into two sets based on their 
apparent symmetry: (1) those for the sky-position and (2) those for the rest of the parameters. The first set displays 
a strong axisymmetry (i.e., symmetry with respect to translations along the (j) axis) that is missing in the other set. 
The reason behind this divide is that the Fisher information matrix and, therefore, the parameter covariance matrix, 
is predominantly block diagonal, with two blocks. The first block comprises just the sky-position (6>, (j)) variance- 
covariance and the second one comprises the variance-covariance of the rest of the parameters. The elements outside 
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FIG. 11: The template- norm- squared for the A and E TDI variables as functions of sky position: I. for the A TDI data 
combination (left), 11. for the E TDI data combination as a function of sky position (center), and III. for the the network 
comprising both A and E TDI data combinations. All plots are presented for the following source parameters: = 7r/4, '0 = 
7r/3, = 3 mHz}. 




FIG. 12: A single arm's template-norm-squared. All plots are presented using the following parameters: {l = 7r/4,'0 = 
7r/3, = 3 mHz}. 

of these two blocks are not all zeros, but are much smaller than the elements within the blocks. Consequently, any 
error in a parameter in one set influences that in another parameter from the same set much more strongly than in 
the other set. This explains why the patterns in each case of the second set are so similar among all parameters within 
that set, as is manifest in Fig. El Also, the reason why the sky-position error plots are predominantly axisymmetric 
while the other ones are not is because while the latter are affected by the beam-pattern functions alone, the former 
are affected by them as well as their parameter derivatives. It is worth noting that the deviation from axisymmetry 
is still apparent in Figs. ^1 however, this effect is significantly reduced with respect to the other parameters. 

There is an alternative way of understanding the three "unphysical" cases (i.e., cases (a)-(c)) studied in Figs. IT^ 
and [131 Starting with the last of these, case (c), note that it is identical to the physical case where the chirping is 
negligible. This is indeed the case when the chirp mass is as low as O.SMq. As seen in Fig. ^ for such a case, filtering 
the data with chirp-less templates causes negligible SNR loss. Thus, the case (c) plots in these figures indicate the 
expected random errors in such a search. 

Case (b) resembles the physical case where the Doppler-phase modulation is negligible, which occurs when the 
source frequency is below a milli-Hertz (barring the slight ruining of the sky-position estimates). Since the chirp rate 
was fixed in this plot, to compensate for lowering the frequency to below 1 mHz one must increase the chirp mass, 
which affects the signal amplitude but not the SNR-normalized error shown here. Thus, case (b) is the expected error 



19 



plot for all sources that have frequencies less than about 1 mHz but chirp masses commensurately larger than 5.3Mq 
(such that mI^^u^^^^ = constant). Comparing the plots for cases (b) and (d) then shows that the SNR- normalized 
error tends to increase at low frequencies, more so at the poles than near the ecliptic. 

Finally, case (a) represents the physical case where u < 1 mHz and A^c is very small, such that both Doppler-phase 
modulation and source chirping are negligible. (Note that the sky-position error is almost constant in the LISA band 
below 1 mHz (cf. Fig. E)).) However as Fig. [HI illustrates, this will not occur in the LISA band. 
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FIG. 13: These are plots of the errors associated with determining the sky position measured in terms of a sohd angle spread 
Afis. The four plots are for the four cases described in Fig. All plots are made for the following source parameters: 
{l = TT/4,tp = 7r/3, V = 3 mHz}. 





FIG. 14: These are plots of the errors associated with determining the dimensionless emission frequency parameter TAi/. The 
four plots are for the four cases described in Fig. El All plots are made for the following source parameters: {l — 7v/4,ip = 
7r/3, = 3 mHz}. 



IX. CONCLUSION 

In this paper we studied the effect of Doppler-phase modulation, frequency evolution, and time-delay interferometry 
on the accuracy of estimating signal parameters of low-mass circular binaries. All results presented here were obtained 
for maximum-likelihood estimates of binary signal parameters in a specific network of LISA's TDI variables, namely, 
the AET variables, and for a one-year integration time. We find that Doppler-phase modulation not only improves 
the sky-position estimation accuracy for z/ ^ 1 mHz, but also degrades it in the range O.lmHz ^ u ^ 1 mHz. It also 
increases the error in the chirp parameter by a factor of 2 and the frequency by a factor of 1.5 at z^ :^ 10 mHz. It 
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FIG. 15: Plots of the errors associated with the remaining parameters: AH/H, T^Az>, and AQl- They all correspond to case 
(d) described in Fig. Both Doppler-phase modulation and frequency evolution are included here. All plots are made for the 
following source parameters: {l = 7r/4,'0 = tt/S^v — 3 mHz}. 

leaves the accuracies of the remaining parameters, namely, the source orientation and the signal amplitude mostly 
unaffected. 

The exclusion of the chirp (or frequency evolution) parameter affects the SNR loss by different amounts for binaries 
with different chirp masses. For instance. Fig. El shows that for Mc ^ 5.3 M©, such an exclusion causes negligible 
effect on the SNR for u ^ 1 mHz. Since in this low-frequency band the estimation accuracy of all signal parameters, 
except u and z>, is affected by its inclusion, from a few percent for the sky-position error to a factor of 6 for the source 
orientation, it is advisable to drop the chirp parameter from the search templates for these sources. On the other 
hand, since the SNR loss from excluding the chirp parameter cannot be ignored when v ^ 3 mHz, for all but the 
lowest chirp masses of interest, one has to accept the parameter-error consequences for these high frequency sources. 
Luckily, a majority of the parameters remain unaffected in this high-frequency band, except u and z>. Thus, if one 
were to turn off or on the inclusion of the chirp parameter judiciously, so as to get the most accurate estimates possible 
for then for a given chirp mass, the error profiles in the frequency plots would be such that they are 

low (i.e., at the levels shown in Fig. 2t) at the low frequencies, but rise at the mid- frequencies, and fall off at the 
high frequencies as in Fig. The biggest casualty arising from searching for the chirp at high frequencies is z>: its 
error can suffer by as much as a factor of 5 at frequencies greater than 3 mHz (see Fig. [SJi). A similar conclusion 
was arrived at earlier by Takahashi and Seto 16] while working with the Michelson variables, although their exact 
numbers were somewhat different. Clearly, this will have an adverse effect on the resolvability of the confusion noise 
above 3 mHz. 

Comparing the frequency plots we obtained here to those found in Ref. [TgI, one striking aspect is the reduced 
levels for the SNR-normalized errors we find for most parameters at all frequencies. Its origin lies in the two principal 
differences between our studies, namely, the choice of the data variables (TDI versus Michelson variables) and the 
approximations made in Ref. p^. For instance, the value of AH/H presented in is its average over LISA's 
band, whereas we derive the values of AH/H for different frequencies. Also, we do not use the long- wavelength 
approximation, which has been shown to have significant effects at frequencies above 3mHz p^ . 

One of the many interesting astrophysical sources that LISA can target is a stellar-mass binary containing a black 
hole. It is unclear how many such objects are present in our galaxy. These objects will likely have chirp masses in 
excess of 2Mq. Figure^ shows that for u ^ 2 mHz, detecting them will require inclusion of the chirp parameter in 
the search templates. This, however, will result in large inaccuracies in all of its parameters in the range 2-3 mHz. 
That in turn will hamper our ability to discern their signals from those of the confusion noise, unless their SNRs are 
sufficiently large. A large SNR can result either from the proximity of a source or a large chirp mass or both. A large 
chirp mass, however, will not always serve this purpose because the 2-3 mHz range actually widens towards lower 
frequencies with increase in chirp mass. This is because the point at which m{'d^'d') starts deviating from unity in 
Fig. shifts to lower frequencies with increasing Aic This identifies a problem with LISA's ability to detect stellar 
mass black hole binaries, about which very little is known in terms of their formation and demographics. 

Note that the above conclusions will change, in general, as the integration time is increased. For example, a ten- year 
integration time will not only improve the accuracy with which the chirp parameter will be determinable, but also 
reduce the adverse effects that its inclusion in search templates has on the estimation accuracy of other parameters 
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Some changes may also appear in the SNR-normahzed frequency and sky plots if one were to make them for a 
different TDI network. However, the major change in that context will arise from the different sensitivities that the 
different TDI networks have to different sky-positions, thus, affecting the absolute (unnormalized) error values. These 
aspects will be explored in more details elsewhere |43j . 
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